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Abstract 

Galactic cosmic rays consist of primary and secondary particles. Primary cosmic 
rays are thought to be energized by first order Fermi acceleration processes at supernova 
shock fronts within our Galaxy. The cosmic rays that eventually reach the Earth from 
this source are mainly protons and atomic nuclei, but also include electrons. Secondary 
cosmic rays are created in collisions of primary particles with the diffuse interstellar gas. 
They are relatively rare but carry important information on the Galactic propagation of 
the primary particles. The secondary component includes a small fraction of antimatter 
particles, positrons and antiprotons. In addition, positrons and antiprotons may also come 
from unusual sources and possibly provide insight into new physics. For instance, the 
annihilation of heavy supersymmetric dark matter particles within the Galactic halo could 
lead to positrons or antiprotons with distinctive energy signatures. With the High-Energy 
Antimatter Telescope (HEAT) balloon-borne instrument, we have measured the abundances 
of positrons and electrons at energies between 1 and 50 GeV. The data suggest that indeed 
a small additional antimatter component may be present that cannot be explained by a 
purely secondary production mechanism. Here we describe the signature of the effect and 
discuss its possible origin. 

PACS codes: 95.30.Cq, 95.85.Ry, 96.40. Cd, 96.40.De 

Keywords: cosmic ray sources, positrons, electrons, supersymmetric dark matter, giant 
molecular clouds 
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1. Introduction 

Over the last 30 years, a number of efforts have been aimed at the study of cosmic-ray 
electrons and positrons with balloon-borne instruments. At energies between about 1 
and 10 GeV, early measurements |l|, found that the "positron fraction" e+/(e+-|-e^) 
is essentially in agreement with a prediction |^ where all positrons are assumed to be of 
secondary origin, and propagate according to the prescription of the simple "leaky-box" 
model of the Galaxy. This is illustrated in Figure la, which shows a compilation of 
measurements []I], 0, P|-P^ of the positron fraction as a function of energy between 0.05 GeV 
and 50 GeV. The leaky-box prediction is shown as a solid curve. At energies around 10 GeV 
and above, as shown in Figure la, several measurements 0, ^ |^ reported a significant 
excess of positrons over the fraction expected from secondary sources. This spurred a 
number of interpretations, ranging from inefficient production of primary electrons at high 



energies [0, |T^ , to hypothetical new sources of positrons [|16| 



The HEAT balloon-borne instrument was designed and optimized to improve the 
accuracy with which cosmic-ray electrons and positrons at energies from about 1 to 50 GeV 
can be detected. The instrument and its performance during two balloon flights in 1994 



and 1995, respectively, are described in detail elsewhere ||13| , |24| - Pq] . A compilation of 
positron fraction measurements is shown in Figure la, where the HEAT results from the two 
flights combined are shown as fllled squares. The overall proton rejection factor achieved 
was nearly 10^. Backgrounds due to atmospheric secondary electrons and positrons were 
estimated by Monte-Carlo techniques, and compared with measured growth curves EA . 



Such backgrounds amounted to 1-2%, and 20-30% of the total electron and positron signals, 
respectively. The uncertainty in the secondary corrections translated in a systematic 
uncertainty of ~ 0.01 in the positron fraction, comparable to the statistical uncertainty; 
however, any systematic error in the correction would affect all data similarly, resulting in 
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an overall normalization shift in the positron fraction distribution, preserving any structure 
observed. 



2. Secondary Production 

The HEAT results shown in Figure la did not confirm the previously-reported rise 
in the positron fraction starting at about 10 GeV. However, the data deviate from the 
predictions of a purely secondary production mechanism in two ways. First, at energies 
below about 5 GeV the positron fraction was in excess of the expectations. For this 



low-energy energy region, another recent measurement [|12[ also reported a positron fraction 
that was significantly higher than measured in the 1960's and 1970's. A possible explanation 
of this effect would come from a solar modulation mechanism that depends on the charge 
sign of the particle and changes from one solar cycle to the next |jlO| . 

The second feature of the HEAT results is an indication of some structure in the 
energy dependence of the positron fraction above 7 GeV. This cannot be easily explained in 
terms of conventional secondary production mechanisms. As shown in Figure lb, a slight 
enhancement in the positron fraction between about 7 and 20 GeV is observed, which may 
suggest a primary source of high-energy positrons. This feature appears in the HEAT data 
from each flight taken separately. Figure lb shows two predictions for interstellar secondary 
production in the energy region of interest. First, the leaky-box prediction is shown as 
a solid red curve. A band of uncertainty in this prediction due to the various uncertainties 
on the parameters of the model and that of the overall normalization is indicated (hatched 
area). In this model, the spectrum of cosmic-ray positrons from secondary sources is 
calculated in the leaky-box approximation from: 
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where {t{E)) is the mean cosmic-ray age at energy E, related to the rigidity-dependent 
mean Galactic escape length, n is the mean density of interstellar nuclei, Pe{E) is the 
rate of production of positrons in interstellar nuclear interactions, and (dE/dt) is the 
rate of energy loss from synchrotron, inverse Compton, bremsstrahlung, and ionization 
processes. The positron fraction is obtained by dividing the predicted positron spectrum 
by the measured all-electron spectrum. A more recent calculation [^, shown as a dashed 
curve in Figure lb, uses a more realistic Galactic diffusion model to predict the positron 
fraction from secondary production. Qualitatively, it predicts the same behavior, a smooth, 
monotonic decrease of the positron fraction without spectral features. The HEAT data 
cannot be well fit by the secondary-production curves of Figure lb. The confidence level for 
the leaky-box prediction is essentially zero (x^=96.5 for 9 degrees of freedom), while that 
for the diffusion prediction is 0.9% (after adjustments to take into account statistical runs 
in the data |2^). Although the band of uncertainty in the predictions is wide, all smooth 
curves within it yield a similarly poor agreement with the data. If the structure seen in 
the data is real, it would indicate the onset of something new, such as an exotic source of 
high-energy positrons. Here we consider several possible models. 

3. Sources of Primary Cosmic-Ray Positrons 

3.1. Annihilating Dark Matter WIMPs 

First, it has been proposed that annihilating Galactic-halo dark- matter WIMPs 
(Weakly-Interacting Massive Particles) are a source of high-energy positrons 



p2|, |23|, |29|, 30|. As most dark matter candidates are Majorana particles, direct 
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annihilation into e"'"e~ pairs is suppressed. In order to account for an observable e"'"e~ line, 
a large total WIMP annihilation cross-section is required. The WIMP density would then 
likely be low and not a major contributor to the present-day cosmological mass density 
||23| . One exception is a model by Kamionkowski and Turner (hereafter referred to as KT) 
20| in which WIMPs with mass m^^ greater than 80 GeV/c^ or 91 GeV/c^ can annihilate 



through resonant production of W"''W^ or Z'^Z° pairs. The resulting electrons and positrons 
are propagated in a leaky-box model. The model predicts enhancements in the positron 
fraction near energies of m^/2 (due to direct decays of the gauge bosons into e^), and 
m^/20 (continuum radiation due to more complex decay chains through intermediate 
production of r^, vr''^, quarks, etc.). If the experimental feature we observe is real, it could 
be a signature for the low-energy continuum radiation peak at around m^/20. Figure 2a 
shows a comparison of the positron fraction we have measured with the HEAT instrument 
and a model prediction including a WIMP-annihilation contribution. On the figure, the 
dashed curve is a baseline secondary production distribution. To this baseline, we add 
a contribution from annihilating dark matter neutralinos pO] of various masses, with an 



amplitude factor left as a fit parameter. This amplitude factor is highly uncertain in the 
model, owing to a combination of large uncertainties in the WIMP annihilation details 
and the astrophysical parameters. The best-fit curve is shown as a solid line in Figure 2a, 
and occurs for a neutralino mass of 380 GeV/c^. The best-fit WIMP source strength is 
1.8 times greater than the estimated amplitude of the effect in the model, well within the 
uncertainties in the prediction. The fit results are summarized in Table 1. The resulting 
confidence level of 74% is markedly better than for fits to strictly secondary production 
models. 

In recent work by Baltz and Edsjo (hereafter referred to as BE) |3^, positron 
production by annihilating dark matter neutralinos is revisited, and a large fraction of the 
Minimal Supersymmetric Standard Model (MSSM) parameter space is sampled. Again, 
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decays and/or hadronization of the annihilation products are simulated, and positron fluxes 
calculated, but a more complex diffusion model than in the KT scenario is used for the 
propagation of the electrons and positrons. Here again, the predicted enhancement in the 
positron fraction is allowed to be renormalized by a factor that is obtained by fitting to 
the HEAT data. Two typical resulting best fit curves are shown in Figure 2a, as dotted 
and dot-dash curves for 336 or 130 GeV/c^ neutralinos, respectively (for details of assumed 



MSSM parameters for these and other models, see |^), and the fit results are summarized 
in Table 1. Once again, an improvement is obtained compared to secondary models, but 
the resulting confidence levels of 22% and 42% are not as high as the best-fit KT model; 
this is mainly a result of the different propagation model used. 



3.2. Pair Creation Near Discrete Sources 

Second, primary positrons could arise when e+e~ pairs are created by electromagnetic 
processes, for instance through the conversion of high-energy 7 rays in the polar cap 



region of Galactic radio pulsars [0. In this model, the positron production rate Pe{E) of 
equation (1) is replaced by: 

Pf^(E) = (l + A;E°-5)Pe(i^), (2) 

where k is given in terms of the Galactic pulsar birth rate 630 (in units of 30 yr), the 
effective time tmax during which the pulsar emits 7 rays (in units of 10"^ yr), the ratio /+ 
of positrons escaping the pulsar per 7 ray produced, and the total interstellar mass M (in 
units of 5 X 10^ solar masses), by: 

^_Q37J0J + Wr_ (3) 

By enhancing the baseline positron fraction from secondary sources with this kind of 
contribution, with k left as a fit parameter, we obtain the best-fit curve shown as a dashed 



line in Figure 2b, with k = 0.15, comparable to reasonable expectations [0. The fit results 
are summarized in Table 1. Although the resulting confidence level of 50% is larger than 
that for purely secondary sources, the shape of the enhancement, a slow monotonic rise 
with energy, is not as compatible with the data as the local-enhancement effect obtained 
with some WIMP-annihilation models. This model predicts that the positron fraction 
should rise with energy beyond 10 GeV, reaching an asymptotic value of 0.5. This could be 
verified with measurements extending to higher energies. 

Another electromagnetic process would be the interaction of very high-energy 7 rays 
with optical and/or UV radiation in the vicinity of discrete sources [^ , resulting in e^ pair 
production. In this scenario, the positron production rate Pe{E) of equation (1) becomes: 

where Et^ = {fn.e(?Y / ^q is the threshold energy for gamma rays interacting with ambient 
photons with characteristic energy eo, x = AE/Eth, T^-y is the optical depth accumulated 
by the gamma ray before escaping the source (a free parameter), and the numerical factor 
is calculated from formulas assuming a gamma-ray power-law index a ~ 2.1. By adding 
a primary positron component from this effect to the baseline from secondary sources, 
assuming various mean values for the parameter eo and allowing the strength of the source 
to remain a free parameter, we obtain the best-fit curve shown as a dotted line in Figure 2b. 
The best fit occurs for eo = 30 eV, as summarized in Table 1, which is in agreement with 
reasonable expectations [0, and requires a relatively weak source strength. The resulting 



confidence level of 75% is once again better than that for purely secondary sources. 
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3.3. Positron Production in Giant Molecular Clouds 

A third possibility is the generation of electrons and positrons in hadronic processes. 
In one model |l^, hadronic cosmic rays can enter and interact within giant molecular 
gas clouds, resulting in the secondary generation of mostly tt^ and K^, which ultimately 
decay into muons, and thereafter into electrons and positrons. Fermi reacceleration due to 
fluctuations in the magnetic field in the turbulent gas could then boost the energy of the e^. 
In this model, if the typical field strength in the cloud is B and the minimum turbulence 
scale is Lmm, a characteristic magnetization momentum p* = eBc/Lmin is defined. Particles 
with momentum greater than p* tend to escape the cloud, so that the spectrum of particles 
accelerated inside the cloud shows an enhancement near p*. For reasonable choices for the 
parameters in the model, it is possible to obtain p* = 10 GeV/c, and a positron fraction 
curve is obtained with an enhancement starting near 10 GeV. If we add to the baseline 
secondary positron fraction such a primary component from giant molecular clouds, and 
allow the strength of the effect to be a free parameter, the resulting best fit is the solid 
curve of Figure 2b. A relatively weak source is sufficient (see Table 1) to fit the data with a 
confidence level of 80%. 

3.4. Other Positron Sources 

Other primary positron sources have been suggested as well. For example, e"'"e~ pair 
production in the magnetosphere of pulsars could be followed by particle acceleration to 
relativistic energies in the pulsar wind driven by low-frequency electromagnetic waves fl^ . 



Or else Z?"*" radionuclei such as ^^Co ejected during a supernova blast, possibly followed by 
shock acceleration in the envelope |^1|, could result in an enhanced high-energy positron 



population. The uncertainties in the models and in the data are such that none of these 
models can yet be ruled out. 
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4. Conclusions 

In conclusion, the HEAT positron measurements indicate a subtle feature, which 
cannot be caused by the atmospheric corrections applied to the data or any other known 
systematic effect. If confirmed by future measurements, this feature could be evidence for 
a new exotic source of positrons, especially at energies beyond about 7 GeV. The exact 
nature of this source cannot be determined until higher-statistics measurements of the 
positron fraction near 10 GeV, as well as extensions of the measurements to energies beyond 
50 GeV, become available. All models proposed and discussed here are unconventional if 
not contrived. Only through detailed studies of the exact shape of the spectral features in 
the positron fraction can we expect to determine which, if any, of the models offers the 
correct description. 
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Table 1: Statistical agreement between HEAT positron fraction results and various primary 
positron source models. 

Model Fit parameter Source amplitude factor Confidence level 



KT WIMPs 


mx 


= 380 GeV/c2 


1.81 ±0.53 


74% 


BE WIMPs 


mx 


= 336 GeV/c2 


11.7 


22% 


BE WIMPs 


^x 


= 130 GeV/c2 


54.6 


42% 


Pulsar 7 rays 




A; = 0.15 




50% 


77 interactions 




eo = 30 eV 


0.0262 ± 0.0076 


75% 


Giant molecular clouds 






0.097 ±0.029 


80% 



The "source amplitude factor" is an arbitrary normalization that indicates the best-fit 
strength of the effect compared to the one predicted by the authors of the model. 
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Fig. 1. — a Compilation of measurements |l|, |^, ^ -13| of the positron fraction between 0.05 
and 50 GeV. The sohd curve is a model calculation ^ assuming that all positrons are from 
secondary sources, and propagate according to a simple Galactic leaky-box model. "HEAT- 



combined" refers to the combination [^ of the data sets from the two HEAT flights, b 
The positron fraction measured with the HEAT instrument, shown on a vertical linear scale. 
The solid curve is a leaky-box secondary model prediction 0, surrounded by an estimated 
band of uncertainty shown as the cross-hatching. The dashed curve is a secondary model 
prediction using Galactic diffusion [^ . 
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Fig. 2. — a The HEAT positron fraction compared with best-fit model predictions with 
an additional positron component arising from annihilating dark matter neutralinos. The 
dashed curve is the baseline solar-modulated leaky-box secondary-production prediction [|10 



renormalized by a factor of 0.85. The solid curve shows an increased positron content due 
to annihilating 380 GeV/c^ neutralinos in the model of Kamionkowski and Turner [^. The 
dotted and dot-dash curves show an increased positron content due to annihilating 336 or 
130 GeV/c^ neutralinos, respectively, in the model of Baltz and Edsjo |30|. b The HEAT 
positron fraction compared with best-fit model predictions from astrophysical sources of 
positrons that are in addition to secondary production mechanisms. The dashed curve is the 
positron enhancement resulting from high-energy 7 rays converting to e"'"e~ pairs near the 



magnetic poles of pulsars [IS]. The dotted curve represents a positron enhancement due to 



high-energy 7 rays interacting with low-energy optical or UV photon fields ||T6[. The solid 
curve shows the enhancement from cosmic-ray interactions within giant molecular clouds 



